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CCD Imaging of the Inner Coma Jets of Comet P/Halley 

James Boswell & David W. Hughes. Department of Physics, University of Sheffield, Sheffield, S3 7RH, 


Abstract. In this paper we analyse the inner coma section of a CCD image of comet 
P/Halley taken at 1807 UT on 1986 March 13th using a C 2 filter (wavelength 5000 to 5200 A 
half maximum) with the 3.8 m Anglo Australian Telescope at Siding Springs, Australia! 
Atmospheric turbulence leads to a spreading of the image detail and this produces a blander 
image of the inner coma region with a slower radial decrease of brightness in comparison to 
the unaffected image We remove this smearing by utilising the point spread function of a 
star on the same CCD image. Jets were then revealed by removing the average background. 
Analysis of the jet structure enabled us to estimate the lower limit of the parent molecule 
velocity. This is found to be 0.3 km s 


Introduction 

JeSun sublimation of ice from active areas of cometary surfaces leads to the ejection 
f h US i!" u 11 ^ velocit y jets. These jets spread out to form a large, quasi-spherical cloud, known 

as the coma, which is, close to the Sun, typically 10 5 km across. The coma surrounds the small solid 
nucleus and is responsible for scattering the majority of the incident sunlight. 

1 his paper investigates the structure of these emissions using CCD images of comet P/Halley taken 
with a C 2 filter (wavelength 5000 to 5200 A, half maximum) at the prime focus of the 3.8 m Anglo 
Australian Telescope at Siding Springs, New South Wales. The image coma region corresponds to a 

SxeT s e i45 to 203 onl “ th * f P .‘ ane ° f ? e . sk f Sf sing throu g h ^e cometary nucleus, encompassing 
pixels 145 to 203 on the horizontal axis and pixels 264 to 322 on the vertical axis. In displaying the data 

ronw°r et T>f S rrn U t0 p ot , b f lghtness 35 a <hei S ht ’ above the two dimensional pixel array. Using this 
extendi _ the , CCD ‘™ ag f e of tbe . comet takes the form of a ‘mountain’ and it becomes convenient to 
extend the analogy and refer to brightness as height and small local intensity peaks as ‘hillocks’. 

Overcoming the problem of atmospheric blurring of the CCD image 

turbulence grouad - ba f d agronomical images is distortion caused by atmospheric 

IhZw u- i? effeC , t J ° f tbe atmo fPhenc turbulence is to independently smear out each of the point 
BaZl lrff would ma ke up the undistorted image, into a twodimensional brightness pattern 
Background stars can be regarded essentially as point sources. Hence, in an ideal image they should 
appear as a single bright pixel. However, in reality, atmospheric blurring causes their image to be spread 
the owr a lT any p | xels f , H f nce > ,n an ir P a ge of the comet they appear as small hillocks on the surface of 
the k if j nta *, r ? tbat re L pre l sen v S the cometary brightness distribution. The height of a hillock above 
the background indicates the brightness of the star that it represents, and in order to ensure energy 

eq!lM7h a e heUhfof to 3 heigbt ° f the hiUock (j.e. the sum of all the hillock DJV-values) rmfst 

the same T lv . alen , t pomt * n ^ he ^distorted image. However, the shape of all the hillocks are 

regardless of , the brightness of the star that they represent. Clearly, every point element in the 
• . t ? rted >7 a ge has been affected in exactly the same manner as the background stars. Thus the final 

image is simply the superposition of all these resultant smeared intensity distributions In mathematical 
to r the ‘STL* d ‘ stnbutl0 . n is convolved with a point spread function (PSF) that corresponds 

Fourier ^ FT*’ 1 the f e t>cally possibk to reconstruct the undistorted image using 

f J Theory to deconvolve the PSF from the CCD image. Unfortunately a secondary 
Transform method is to reduce the signal-to-noise ratio of the image. In the case 
fundfstort H H ^ ‘TP U WaS f ° U ? d ^ hat ’ the S / N ratio was reduced to such an extent that the signal 

louffht ?n beCai ¥ C ?P p etely swamped b y the noise - Thus, an alternative approach must be 

sought in order to deconvolve the comet image successfully. 

Image Deconvolution 

The image of the comet can be more easily processed if the underlying Cartesian pixel array is replaced 
by an equivalent polar array which has the ‘summit’ of the mountain at the orig^7 an77radiaf pixel 
shonidV ha J 1S th f SamC “ , the ° u nglnal Cartesian spacing. The azimuthal spacing of the polar array 

the DW ' vat of alfofri! h i ’? T ' OSS ° f re f oIution at the edges of the new imaged veraging 
e UN values of all of the polar pixels for a complete azimuthal sweep then reduces the imatre to a 

nr^fitT enS1 |° na brightness profile. Because of the inherent circular symmetry of the image, this average 
profile is a close approximation to any radial slice that is taken through the mountain. The same procedure 
may be applied to produce an brightness profile for a background star within the same field of view as 
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the comet. An approximate fit of the brightness profile of the comet may be obtained from the equation 
DN(R\ — R~ n (R is the radial pixel number, where the central pixel in the image is num , ( ) 

UhLs »d W positive real number). However, n varies with Rover .narrow '“g 

S 1 "ell profile" havin/values in the range 0 to 2. We can therefore . generate i a. 
brightness profiles for the cometary image by varying n between these limits. These hypotheti p 
with the brightness profile generated for the background 
The percentage change in brightness found on convolving the hypothetical profiles with the Fbt y 

theii be calculated^ ^ g ig < 1% for pixels whe re R > 10 for all values of n, thus pixels with 

p > 10 ma y be ignored in the subsequent selection of the value of n that will ultimately be used. For 

fllpixels the change in brightness is negligible for 0.25 < » ■ < 0 40 - 

prt a value of n which produces a brightness profile which best matches that ot pixels wim ^ 
Thfv^ue of n which provides the best fit for the data is 0.37, so a hypothetical Witness profile was 
j - i . * d at ( d— 0.37 this Drofile is then convolved with the PSF. We may 

L7*dS?v Tp «nn between ?h. original hypothetical profile with , „ = = <U7 
and the mw profile after convolution. Each pixel of the original comet image can then be divided by the 

Sri|htisSs a wLe Ca found to^ee to'wi^ri^oveT Tmale^The deconvolved image is shown 
in figure 1 . 



Figure 1. The deconvolved CCD image. The base grid represents the pixel array 
of the CCD image of comet P/Halley. The pixel numbers are shown. The vertica 
axis represents the deconvolved pixel DN value. Note the sharp peak at the central 
position i.e. at the nucleus. The base is at a DN ‘height of 2460 with the peak at 

12220 . 
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Fig ure 2 . The deconvolved CCD pixel array with the average background removed, 
this is done by converting the pixel array from cartesian to polar coordinates, then 
finding the mean DN value over all d for each value of R, and subtracting this from 
all pixels at distance R. Note the two raised regions (jets). The C 2 enhancement 
corresponds to, at maximum 9%. 


Enhancing Image Detail 

?Wv2 e (, imag u ha t be€I [ deconvolved, the average background brightness was removed This was 
polar c^r/inTtes^aftei^ d^onv^^fon^fro^^h^p^xeF'accor^in^^ t0 

imase (figL 2) - ™ s 

Calculating Jet Velocity 

J he structure ofl the jets it is possible to make an estimate of the outflow velocity assuming 

Wrnmmmmmi 
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Unfortunately the rotation period of P/Halley is still the subject of investigation. Itoh (1985) first 
discovered a 2 2 day periodicity, but within a year Millis and Schleicher (1986) added a 7.4 day periodicity. 
Subsequent analysis tupports^both periods, leaving the issue unresolved. For the purpose of calculating 
the outflow velocity of the jets, ti Gl the 7.4 day period provides a bwer limit va km • If the 

2.2 day period is proved correct, however, this limit may be raised to v G > d.d km s . 



Figure 3. Both figures are essentially figure 2 redrawn with the abscissa representing 
the azimuthal distance 9, and the ordinate representing the radial pixel distance R. 
The grey scale used in the upper plot is the same as that used in figure 2. The lower 
plot has contours at -450, -250, -50, 150, 250 and 350. The degree of curvature of the 
jets can be assessed from their gradient on this plot. 
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